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Abstract 

I The evolution of a universe modelled as a mixture of generalised Chaply- 

' gin gas and ordinary matter field is studied for a Robertson Walker type of 

CN . spacetime. This model could interpolate periods of a radiation dominated, 

matter dominated and a cosmological constant dominated universe. Depend- 
I ing on the arbitrary constants appearing in our theory the instant of flip 

changes. Interestingly we also get a bouncing model when the signature of 
one of the constants changes. The velocity of sound may become imaginary 
under certain situations pointing to a perturbative state and consequently 
the possibility of structure formation. We also discuss the whole situation in 
5^ . the backdrop of wellknown Raychaudhury equation and a comparison is made 

W)| with the previous results. 



O 



X 



PACS : 04.20, 04.50 +h 



! KEYWORDS : cosmology; accelerating universe; chaplygin gas 

> 
00 

\o 
cn 

in : 1 Introduction 
o 

Three discoveries in the last century have radically changed our understanding 
of the universe - as opposed to the idea of Einstein's static universe Hubble and 
Slipher(1927) showed that it is expanding. Secondly CMBR as also primordial nu- 
cleosynthesis analysis in the sixties point to an initial hot dense state of the universe, 
which has been expanding for the last 13.5 Gyr. Finally, if we put faith in Einstein's 
theory and FRW type of model then as standardised candles type la supernova sug- 
gest [1] that the universe is undergoing accelerated expansion with baryonic matter 
contributing only five percent of the total energy budget. Later data from CMBR 
probes [2] also point to the same finding. This has naturally led a vast chunk of 
cosmology community to embark on a quest to attempt to explain the cause of the 
apparent acceleration. The vexed question in this field is the possible identification 
of the processes likely to be responsible for triggering the late inflation. Researchers 
are plainly divided into two broad groups - either modification of the original Ein- 
stein's theory or introduction of any exotic type of fluid like a cosmological constant 
or a quintessential type of scalar field. But the popular explanation with the help of 
a cosmological constant is beset with serious theoretical problems because absence 



^ Sree Chaitanya College, Habra 743268, India and also Relativity and Cosmology Research 
Centre, Jadavpur University, Kolkata - 700032, India , e-mail: dibycndupanigrahi@yahoo.co.in 

^ IGNOU Convergence Centre, New Alipore College, Kolkata - 700053, India and also Rel- 
ativity and Cosmology Research Centre, Jadavpur University, Kolkata - 700032, India, e-mail : 
chat ^uj it 1 @yahoo .com 
Correspondence to : S. Chatterjee 



of acceleration at redshifts z > 1 implies that the required value of the cosmological 
constant is approximately 120 orders of magnitude smaller than its natural value 
in terms of Planck scale [3]. As for the alternative quintessential field [4] we do 
not in fact have a theory that would explain, not to mention predict, the existence 
of a scalar field fitting the bill without violating the realistic energy conditions. 
Moreover we can not generate this type of a scalar field from any basic principles 
of physics. Other alternatives include k-essence [5], tachyon [6], phantom [7] and 
quintom [8]. So there has been a resurgence of interests among relativists, field 
theorists, astrophysicists and people doing astroparticle physics both at theoretical 
and experimental levels to address the problems emanating from the recent extra 
galactic observations without involving any mysterious form of scalar field by hand 
but looking for alternative approaches based on sound physical principles. Alterna- 
tives include, among others, higher curvature theory, axionic field and also Brans- 
Dicke field. Some people attempted to look into the problem from a purely geomet- 
ric point of view - an approach more in line with Einstein's spirits. For example, 
Wanas [9] introduced torsion while Neupane [10, 11] modifies the spacetime with 
a warped factor in 5D spacetime in a brane like cosmology and finally addition of 
extra spatial dimensions in physics as an offshoot of prediction from the string the- 
ory [12, 13, 14, 15, 16]. While this torsion inspired inflation has certain desirable 
features the problem with Wanas' model is that the geometry is no longer Rieman- 
nian. Further a good number of people [17, 18, 19, 20] have done away with the 
concept of homogeneity itself and have argued that accelerating model and conse- 
quent introduction of exotic matter field have to be invoked only in FRW type of 
cosmology. In Tolman Bondi like inhomogeneous model the apparent dimming of 
the signals may be explained as a consequence of inhomogeneous distribution of 
matter. 

While the above mentioned alternatives to explain away the observed acceleration 
of the current phase have both positive and negative aspects the one that caught 
the attention of a large number of workers is the introduction of a Chaplygin type 
of gas as new matter field to simulate a sort of dark energy. The form of the 
matter field is later generalised through the addition of an arbitrary constant as 
exponent over the mass density and is generally referred to as generalised chaplygin 
gas(GCG) [21, 22]. Though it suffers from the serious disqualification that it vio- 
lates the time honoured principle of energy conditions its theoretical conclusions are 
found to be in broad agreement with the observational results coming out of gravita- 
tional lensing or recent CMBR and SNe data in varied cosmic probes [23, 24]. This 
is generally achieved through a careful maneuvering of the value of the newly intro- 
duced arbitrary constant. To further fine tune the match between the theory and 
the very recent observational fallouts the GCG is again modified via the addition 
of an ordinary matter field, which is termed in the literature as modified chaplygin 
gas(MCG) [25, 26]. The viability of such scenarios has been tested by a number of 
cosmological probes, including SNe la data [23, 24], lensing statistics [27, 28, 29], 
age-reshift tests[30], CMB measurements [31], measurements of X Ray luminosity of 
galaxy clusters [32], statefinder parameters [33]. In our previous work [31] we have 
studied Chaplygin gas model in inhomogeneous space time. In the present work we 
have revisited the dynamics of the FRW model taking MCG as matter field and 



have tried to discuss some as yet unexplored region and have got some interesting 
results. We have organised the paper as follows: In section 2 the mathematical 
formulation is given and we have ended up with a hypergeometric solution and also 
an effective equation of state as p = yV(t)p in section 3. So depending on initial 
conditions our model mimics both AC DM and quiessence models and the evolution 
is also shown graphically. We have also made some detailed discussion on acoustic 
wave in our model and find that all possibilities like less/greater than light velocity 
and even imaginary values exist in our model. Relevant to mention that imaginary 
sound velocity is not that much discouraging in this context because it gives rise to 
perturbation and consequent structure formation [35]. The interesting thing in our 
analysis is that we have taken the first order approximation of the field equation 
as key equation and subsequently found out the exact solutions. We are not aware 
of attempts of similar kind in the past literature. Moreover it is also found that 
if an arbitrary constant appearing in our solution be taken negative the cosmology 
bounces back from a minimum. We have also made a detailed analysis of flip time 
both analytically and graphically in this section. In section 4 these conclusions are 
checked in the framework of well known Raychaudhury equation. The paper ends 
with a discussion in section 5. 



2 Field Equations 

We consider a spherically symmetric homogeneous spacetime given by 

ds^ = df -a\t) {dr^ + r^dn^) (1) 
where the scale factor, a{t) depends on time only. 

A CO moving coordinate system is taken such that = 1, = (z = 1, 2, 3) and 
g^^u^Uy = 1 where Ui is the 4- velocity. The energy momentum tensor for a dust 
distribution in the above defined coordinates is given by 

Ti: = {p + p)6X - pS'u (2) 

where p{t) is the matter density and p{t) the isotropic pressure. 
The independent field equations for the metric (1) and the energy momentum 
tensor (2) are given by 

3^ = P (3) 

d 

2- + - = -p (4) 
a 

From the the Bianchi identity we get for the homogeneous model the conservation 

law 



V^T^'^ = 

which, in turn, yields 



(5) 



p + 3-(p + |>) = (6) 
a 

At this stage we assume that we are deahng with a Modified Chaplygin type of 
gas (MCG) obeying an equation of state 

V = Ap-^ (7) 

where A, B and a are constants. The exponent a, from most observational con- 
straints, hovers around unity [36] and the constant A ranges from 1/3 to zero. 
Similarly the positive definite constant B is also not exactly arbitrary. In the equa- 
tion (7) when the last two terms start to be of the same order of magnitude the 
pressure vanishes. In this case the fiuid has pressureless densitypo corresponding to 
some cosmological scale Oq given by po = P"'''^(oo) = ^- Many variants of Chaply- 
gin gas model have come up in the literature and the equation (7) refers to what 
is generally known as the Modified Chaplygin gas model (MCG) [25, 26] such that 
A = gives generalized model (GCG) [31] and if in addition a = 1 one recovers the 
original model. Moreover, the first term on the r.h.s. of the equation (7) gives an 
ordinary fiuid obeying a barotropic equation of state (EoS) so that we here are es- 
sentially dealing with a two fiuid model. Further, the equation (7) points to an EoS 
that interpolates between standard fiuids at high energy densities and Chaplygin 
gas fluids at low energy densities. In the 4D framework the dynamics of the MCG 
model has been studied in the reference [25, 26] and a perturbative study looking 
for some generic features is carried out in [37]. On the other hand Fabris et al [38], 
in an interesting work, have used a perturbative analysis to confront observational 
data within this model and taking the particular case of power spectrum observa- 
tional data have concluded that the recent data restricts the value oi A < 10"^ 
such that the GCG is recovered and the MCG is almost ruled out. Moreover, in 
the case of MCG model recent supernova data seem to favour negative values of 
the parameter a [39]. When one attempts to address issues concerning structure 
formation the study of cases with negative values of a becomes more sensible since 
this implies imaginary sound velocities, hence plagued with the possibility of in- 
stabilities [40, 23, 24]. On the other hand it has been argued that a > 1 is also 
plausible [41]. With the help of equations (6) & (7) a little mathematics shows that 
an expression for density comes out to 



3(1 + A){l + a) [ _^a3(i+^)(i+")-i da + c 
J 1 + A 
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l + a 



where c is an integration constant. The above equation (8) yields a first integral 

as 
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l + a 



Plugging in the expression of p from equations (3) and (9) we finally get 



(9) 
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(10) 
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(a) A = 0.5 
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Figure 1: The variation of a{t) and t for different values of a. The graphs clearly 
show that the scale factor a(t) increases in greater rate for smaller a. 



Using equation (10) we have drawn the figure-1, where the evolution of a{t) with 
t is shown. This figure shows that as a increases the rate of change of scale factor 
decreases. A cursory look at the equation also points to this type of variation of the 
curves. In fact the equation (10) suggests that with a, ^ becomes fiatter. 

At this stage if we consider A = (figure-lb) we get from the equation (7) 
that as a increases (0 < a < 1) the first term on the r.h.s. reduces to zero while 
the second term decreases for a particular p. We get identical results for A = 
3(1+") (figure-lc). Now, p being a small fraction at the late stage of evolution of the 
universe any increase in its value in the exponent finally increases its magnitude so 
that the pressure becomes more negative, which, in turn drives the expansion more 
vigorously. 

Our analysis is based on different sets of observational data. By using a large 
sample of milli-arc second radio sources recently updated and extended by Gurvits 
et al [12] along with the latest SNel data as given by Reiss et al [13], Alcaniz and 
Lima [11] showed that the best fit data for these observations are Bs = 0.84 & 
a = 1.0 (UDME) and 5, = 0.99 & a = 1.0 (CGCDM), where = ^, where 
Pq is the present density of the Chaplygin gas. In another work Lima et al [36] 
showed at 95% confidence level by the BAO (Baryon acoustic Oscillation) and Gold 
sample analysis, the range of a is 0.9 < a < 1 while BAO & SNLS analysis provides 
0.94 < a < 1. Both the results predict a to be nearly equal to unity. In contrast 



to this result Fabris et al [3(S] as pointed out earlier ruled out the existence of A 
for the MCG model in the context of power spectrum observational data. In this 
context the relations A = and A = „/~°^ seem interesting. When the value of 

l+a 3(l+a) 

a is nearly equal to unity there remains a tiny value of A, which is not exactly in 
line with the work of Fabris et al [38]. Lastly Lu et al [41] gives for the MCG best 
fit data A = -0.085 and a = 1.724 in the light of 3 yr WMAP and SDSS data. 



3 Cosmological dynamics 

It is very difficult to get the exact temporal behaviour of the scale factor, a{t) 
from the equation (10) in a closed form because integration yields elliptical solution 
only. However, the equation (10) does give significant information under extremal 
conditions as briefiy discussed below. 



Deceleration Parameter: 



At the early stage of the cosmological evolution when the scale factor a{t) is 
relatively small the second term of the last equation (10) dominates which has been 
already discussed in the literature [25, 26]. So we will be very brief on this point. 
From the expression of the deceleration parameter, q we get 



l^i + ^^ 
2 2p 



q = - — - = — [n ^)-i = - + -- (li; 

where H is the Hubble constant. With the help of the EoS given by (7) we find 

1 + 3y4 35 1 



q - 

which via equation (9) gives 
1 + 3A 
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Figure 2: The variations of q and p for different values of a. Here A = 0.001. 



As the universe expands p decreases with time such that the second term in the 
equation(12) increases pointing to the occurrence of a flip when the density attains 

a critical value given by p = pfup = (j^^) This flip density pjup depends on 
the exponent a such that at the larger value of a the flip density decreases, i.e., flip 
occurs at a lower density, i.e., it occurs at a later time. We know from Lima's [36] 
that the value of a is restricted to 0.9 < a < 1 such that acceleration is a recent 
phenomenon. This result is encouraging. As discussed in the last Section Lu et al 
[41] argued that a > 1 also conforms to the observational analysis. This finding is 
particularly relevant to our case in the sense that higher values of a signify a lower 
Pfiip, i.e., more recent accelerating phase. The above analysis in conformity with 
the nature of g ~ p curve in figure-2. 

CASE A : At the early stage when the scale factor, a{t) is very small the 
equation (13) reduces to 

Q = (14) 

Evidently the deceleration parameter has contribution from the baryonic matter 
content only such that, q mimes the ordinary fluid behaviour with magnitudes 1 
and I for radiation and dust respectively as in a FRW model. When A = —1 the 
equation (14) gives, q = —1 evolving as a AC DM model. 

CASE B : In earlier works [25, 26] authors utilized the above equations to 
find an equation of state at the late stage of evolution as, p = {a + {1 + a)A}p. 

Using the equations (7) & (9) straightforward calculations yield an effective EoS 
at the late stage of evolution as 



p = p[A- 5p(i+")] 

where 



{1 + Af 



B ^3(l+A)(l+a) 



p = W(t)p (15) 



which is a function of time only. This is clearly at variance with the earlier works 
of [25, 26] where the effective EoS shows no time dependence. We also find that 
at the late stage of evolution as a{t) — )■ oo, W(t) — ?■ —1 so we asymptotically get 
p = —p from this Chaplygin type of gas, which corresponds to an empty universe 
with cosmological constant such that the equation (11) implies that the deceleration 
parameter, q reduces to —1. Interestingly W{t) always remains greater than —1, 
thus avoiding the undesirable feature of big rip. In this context we call attention 
to a recent work of Z. K. Guo and Y. Z. Zhang [45] where a new variant of CG is 
taken in the form of 

P--^ (17) 

P 

where unlike the original CG, B is taken as a function of the scale factor a{t). For 
mathematical simplicity they assumed B{a) = Bqu""" where Bq and n are constants 



and n < A and Bq > 0. They finally end up with a constant equation of state 
parameter 



w = -l+n/6 (18) 

We find that n = corresponds to the original Chaplygin gas model which inter- 
polates between a universe dominated by a dust and DeSitter era. Moreover n > 
corresponds to a quiessence dominated and n < to a phantom dominated model. 
In our case we, however, get here a time dependent equation of state parameter 
which always avoids the undesirable phantom like behaviour. 



To end up a final remark may be in order. In an earlier work the present authors 
[12, 13] in the framework of {d + 4) homogeneous spacetime studied the scenario 
with an EoS given by equation(7) but generalised to extra dimension. Using an 
ansatz b{t) = a{t)~^ where a{t) and h{t) are 3D and extra dimensional scale factors 
and m is a constant has led us, at the late stage, to an EoS p = wp. The expression 
for the w is found to be 

2dm{'m + 1) 



w 



k 



(19) 



where k = dm?{d — 1) + 6(1 — dm), is a constant. Unlike the usual 4D cases 
(see for example [16]), here w ^ —1. Obviously this is due to the presence of extra 
dimensions in the above relation. In 4D case (d = 0) w; = — 1 and a AC DM model 
is the only possibility. In general the magnitude of w is parameter dependent and 
presents varied possibilities. When m = 0, i.e. a{t) is a constant we again get back 
the 4D case. When m > 0, w < —1; So a phantom like cosmology results with the 
occurrence of 'big rip' etc. But the cosmology becomes physically interesting when 
— 1 < m < such that > w > —1 and we get a quiessence type of model [47, 48]. 

The variation of W(t) with the scale factor a{t) for different values of a are shown 
in the figure-3. We have considered three cases: the constant value of A = 0.5 is 
chosen for the figure 3a, on the other hand we have chosen the relation A = and 
^ = 3(\"!^^) for the figure 3b and 3c respectively. All the graphs clearly show that 
the scale factor a{t) increases as W(t) becomes more and more negative. 

Since, —1 < W(t) < 0, the relation shows that W(t) can never be less than — 1, 
a good sign. Otherwise there will be a phantom stage. In quintessence model W(t) 
starts from zero and then reduces to —1. 



Acoustic wave : 



In this case the expression of the velocity of sound Vs with the help of equation 
(15) will be 

= I = -4(1 + a) - = A + .(1 + A) |i - ^;3a.t)al) } (20) 



Using equations (11) and (20) we get the expression of the deceleration param- 
eters. 




Figure 3: The variations of Wit) and a(t) for different values of a. All the graphs 
clearly show that the scale factor a(t) increases as W(t) becomes more and more 
negative. 



,= l + ?P = i + M(i±^/l_^^^ (21) 



2 2p2 2a \ A{l + a 

We have considered three relations for a as a = 1, a = and a = to 
study the above equations in a more transparent manner. From the observational 
point of view it is seen that the value of a is nearly equal to the unity. As pointed 
out earlier Fabris et al [3S] studied and ruled out the constant A. However, our 
investigations differ and in a sense more general than Fabris et al [38] in that we 
have allowed a small value of A for 0.9 < a < 1 [36]. When a > 1, we get the 
negative value of A which also is in agreement with some observational result [11]. 

I. {a = 1): 

From equation (20) we get 

''J = '4+(l + ''){l-l^} (22) 

The equation (22) shows that for A = 0, Vs is always less than the velocity of light 
Vc and can not be imaginary. For any other values of A (the limit oi Ais < A < 1) 
the velocity of sound Vg may not be less than Vc- Now with the help of the equation 
(20) we calculate the condition for Vs < Vc which is 



On the other hand A there may be a possibihty for Vs > Vc, but our discussion 
is restricted only to the late stage of evolution where the scale factor a(t) is large 
enough. In this context we get Vg > fc at the late stage of evolution. The above 
phenomenon is shown graphically in the figure-4a. 
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Figure 4: The variation of Vs with a{t) for different values of A h a. In the fig.- (a) 
Vg ^ Vc only for A =0, hut in the fig. -(h) Vg ^ Vc for any value of a (since Q < a < 1) 
or A (since Q < A < 1). In the fig-(c), for a < 1 gives Vg < Vc and for a > 1 shows 

Vg > Vc. 



II: {a = '■ From equation (20) we get 

vt = 1- 



1 + aYBa^ 



(24) 



It is evident from the equation (24) as well as from the figure-4b that the velocity 
of sound Vg always less than the velocity of light Vc for any value of a ( since 
< a < 1). For a = 1, A = and we get back the situation depicted in the 
figure-4a, however, for any other values of a (0 < a < 1), A 7^ 0. Some observations 
predict that the value of a is nearly equal to 1 [3(i]. So for the maximum permissible 
value of a, Vg should be always less than the velocity of the light Vc- But in the 
previous case ( for equation (22)), there may be a possibility that the Vg is greater 
than Vc [19] for the high value of a{t) i.e., at the very late stage of evolution. We 
also get the same conclusion from the equation (23). 
Ill: ia='^): 



From equation (20) we get 



a 



2a{2 



a 



3(1 + «) 3(1 + a) 



1 - 



2c(2 + a) 



35(1 + a) a6(2+") 



(25) 



For late universe when a{t) is very large, the above equation reduces to ~ 
When a = i.e., we get back to the AC DM model, in this case = | and for 
a = 1 (in this case A = 0), i.e., for pure Chaplygin gas model, Vg = l which imply 
that for < a < 1, Vg < Cg. Again for a > 1, we have seen that Vg > Cs, however, 
violates the causality condition. The figure-4c gives similar conclusion. 

Now we discuss the whole analysis of Vg in the context of deceleration parameter 
q using equation (21). To get accelerating universe q should be negative. So the 
condition for accelerating universe involving Wj, is f ^ > ^ + A{1 + a). 
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Figure 5: The variation of q with Vg for different values of A ^ a. 

From the figure-5a it is seen that flip occurs for A = when Vg < Vc- But for 
other values of A, at the time of flip, Vg > Vc- In the figure-5b & 5c at the time of 
flip, Vg is always less than Vc for different values of a or A 

Now our analysis will be restricted within the accelerating phase, i.e., after flip. 
For a = 1, to get acceleration > | + 2A. If we consider A = 0, i.e., when only 
the original chaplygin type fluid is present, > |. So the velocity of sound may or 
may not be greater than the velocity of light. For A ^ 0, the above expression for 
the velocity of sound further implies that A < | in order that Vg < Vc- In this case 
we restrict the limit of A as < A < |. Again when a = to get acceleration, 

> 1 — For a = 1, Vg > | exactly similar to the situation discussed earlier. 

When a = > . For a = or 1, > but for a > 1, Vg > negative 

value. 




CASE C : Distinctly new models unfold itself when we take the arbitrary 
integration constant as c < 0. Here the energy density increases with the scale factor 
mimicing a phantom dark energy model and finally ending up as a cosmological 



constant. We get from (15) that for the matter field to be well behaved the condition 



3(i+A)(i+«) > C{1 + A) 

B{l + a) ^ ' 

need to be satisfied. So a minimal value of the scale factor given by 

'C{l + A) 



a{t)r 



B(l + a] 



1 



3(l+A){l+c«) 

(27) 



This naturally points to a bouncing cosmology at early times. In the past Setare 
[50] analysed these possibilities in a series of work. To sum up we see that the Chap- 
lygin model interpolates between a dust at small a and a cosmological constant at 
large a but this well formulated quartessence idea breaks when a negative value of 
the arbitrary constant is taken. Following Barrow [ 1 6] if we reformulate the dynam- 
ics with a scalar field ( and a potential V to simulate the Chaplygin cosmology, 
we find that a negative value of c implies that we transform ( = i'^. In this case 
the expressions for the energy density and pressure corresponding to the scalar field 
show that it represents a phantom field. 
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CASE D : 

As we are considering a late evolution of our model the last term in the equa- 
tion(lO) is almost negligible compared to the second term and so the findings coming 
from a first order approximation of the equation(lO) may be of relevance. Here we 
find an exact solution of the first order approximation of the equation(lO). Authors 
of this work are not aware of attempts of similar kind in any earlier work. So this is 
clearly a new result. Now from equation (10) we get, as first order approximation 
the equation at the late stage of evolution 

3^=(t^) ...n + ^n^.. (28) 

For economy of space we skip the intermediate steps and write the final solution 

as, 



/ C 1 + A\ 3(l + A)(l + a) 2 f V^S l + 2a 1 , 

a{t) = ( — ^ ^ j sinh3(i+^)(i+") J _(i + yi) 2(i+a) (i _|_ ck)52(i+«) |> t (29) 

From figure-6 we have seen that the role of the parameters (A and a) are just 
opposite to what we observe in figure- 1. A plausible explanation may be the fact 
that unlike the first case only first order terms are present here. So the higher order 
terms in the first case drastically change the scenario and makes their presence felt 
in changing the nature of the curves. At this stage correspondence to our earlier 
works [12, 13] may be of relevance. We have shown that one may get similar form of 
solution in a higher dimensional spacetime if a particular ansatz on the expression 
of deceleration parameter is taken apriori . But the essential difference between the 





Figure 6: The variation of a(t) with t for different values of A ^ a are shown in 
this figure. The graphs clearly show that the rate of increasing of the scale factor 
a(t) increases for greater value of a or for smaller value of A. 



two lies in the fact that while in the earlier work the hyperbolic solution results from 
a particular form of the deceleration parameter here we have to invoke a Chaplygin 
type of gas to get similar comological evolution. 

Deceleration Parameter : 

The equation (29) can be reduced in the following form 

a{t) = aosinh"a;t (30) 



where, ao = {^}^'""'"""^ (1 + A)W^^ , n = 3^,^^ and = ^(1 + 

l + 2a 1 
A) 2(1+-) (1 + Q;)_B2{l+a) 

such that we get from equation (30) 

la 1 — n cosh^wt 

^ = = T:2—— (31) 

H"^ a n cosh ut 

and 

t, = -cosh-if^) (32) 
u \\/n 



showing that the exponent n critically determines the evolution of q. A little 
inspection shows that (i) n > 1 gives always acceleration, (ii) < n < 1 gives the 
desirable feature of flip, although it is not obvious from our analysis at what value 
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Figure 7: The variation of q with t for different values of n are shown in this figure. 
This figures show that (i) n> 1 gives always acceleration and (ii) n < 1 gives flip. 



of redshift this flip occurs. Figure-7 gives the similar conclusion that late flip occurs 
at lower value of n. From equation (31) it further follows that for physically realistic 
values of A and a as positive definite < n < 1 and a fiip is a distinct possibility, 
again it follows from equation (31) and also from figure-8 that the early fiip occurs 
at higher values of n as well as u. Again n and u depend on A and a. It can 
be said that for constant values of u and A, late fiip occurs at higher values of a, 
which have some observational implications that the value of a should nearly equal 
to unity (0.9 < a < 1) or greater than unity. If we observe the expressions of n and 
w, we can not say clearly what values of a and A give the early fiip. But we can say 
about the time of fiip if we consider some special value of a. 



I 2 





0.0 0.2 0.4 0.6 0.8 1.0 1.2 



Figure 8: The variation of tc with n for different values of oo are shown in this figure. 
This figure shows that the early fiip occurs at higher values of n as well as uj. 



Now we consider some special cases. 
I. a = 1: 

We get from equations (31) & (32) for a = 1 we get the expression for deceleration 
parameter 





3(1 + A)- cosh^ (^3(1 + A)lBh} 

q = ^ ^ ^ (33) 

cosh'^lV3il + A)iBn\ 

And the flip time can be calculated from the above equation as 

V3(l+A)454 1^ J 

Ify4 = —1, q = —1, we get the evolution dominated by A with no contribution from 
Chaplygin gas. This also follows from the equation (34) because here tc oo. Thus 
there is no flip as in the de-Sitter model. 

Again using equation (31) for a = we get the expression of the deceleration 
parameter 



1+A 



cosh^ y3(^) ' t 



cosh 



B{l+a) 1 2{l+a) 



From the equation (32) the flip time tc becomes 



(35) 



l+A 1 

tc = ^ (^-^] * cosh-^ (VS) = i^TT^ r I cosh"^ (V3 

V3\ B ) \ ) ^\B{l + a)) V , 



III.( 



a 



1-3A 



1+3A; 

From equation (31) we get 



1 cosh^ \ 

^ 3(1+A) ^"^'^ 1 




1+3.4 ~1 


cosh^ < 




+ A)t(--)^B-r] 





And from equation (32), the flip time tc becomes 

1 + 3A 



^cosh-N/^ii±^ 

v^(l + A)i(i-^)i?^ V 1 + 3A 



(36) 



(37) 



(3^ 



In figure-9, we see that flip depends on a. From the graph it is seen that the 
flip time increase with lower value of a. a{t) — oo, t — > oo in agreement with our 
graph. 



4 Raychaudhuri Equation 

It may not be out of place to address and compare the situation discussed in the 
last section with the help of the well known Ray Chaudhuri equation [51], which in 
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(a) a = 1 (b) a = 

Figure 9: The variation of q with t for different values of A ^ a. The graphs clearly 
show that the flip time is greater for smaller value of A for the fig.- (a) and increases 
for greater value of A or smaller value of a for the fig. -(b). 

general holds for any cosmological solution based on Einstein's gravitational field 
equations. With matter field expressed in terms of mass density and pressure Ray 
Chaudhury equation reduces to a compact form as 

^ = -2(a^-a;^)- V-^(p + 3p) (39) 

in a CO moving reference frame. Here p is the isotropic pressure and p is the energy 
density from varied sources. Moreover other quantities are defined with the help of 

a unit vector v'^ as under 



the expansion scalar 6 = v^;i (40a) 
^2 = aija'^ (40b) 

1 1 1 

the shear tensor (Xij = -{vi-j + vj-i) - -{viVj + VjVi) - -v"^{gij - ViVj) (40c) 

the vorticity tensor Uij = 2 ("^^y' ~ "^i;*) ~ Iji'^i'^j ~ '^j'^i) (40d) 
We can calculate an expression for effective deceleration parameter as 

H + H'^ e 

^ = ^ = -1-3^ (41) 

which allows us to write, 

e^q = 6a^ + 127iG{p + 3p) (42) 

In our case as we are dealing with an isotropic rotation free spacetime both the 
shear and vorticity scalars vanish. 

With the help of the equations (7), (9) & (42) we finally get, 



B 



l + A 



2B c(l + 3A + 3a + 3Aa) 



l + A 1 + a a3(i+a)(i+A) 



(43) 



In original Chaplygin gas where a = 1,A = we get from the equation (44) , 
_ 2nG (^—B + This is exactly similar to what we have found in our earlier 
work [■) I] (vide equation 3.11), when dealing with an inhomogeneous LTB model. 
Now we consider some special cases. 



I. (a = 1): 

From Equation (44) we get 
e'^q = 12ttG 



B 



A 



2B c(2 + 3A) 

a6(i+A) 



A 



(44) 



In this case flip occurs when g = 0, at that time the scale factor a{t) will be 

0,(tflip 



6(1+A) 



(45) 



Now, g < at a{t) > {^{l + A){2 + 3^4)}'*'^+'*) such that acceleration takes place 
in this case. 



II. ( 



a 



l~A 
l+A 



Again using the equation (44) we get 



-B(l + a) + 



4c 



(46) 



Here, at the flip time q = and at that time the scale factor a{t) will 



ait f Hi 



vp) 



Ac 



B il + a) 



B 



(47) 



The acceleration takes place when g < i.e. a{t) > [■§{A + 1)^] 



iii.( 



a 



1-3A 



1+3A, 

From equation (44) we get. 



{ 2(2 + a) 



35(1 + a) ^ 2c 



(2 + a) a2(2+°) 



(48) 



Here, at the flip time q = and at that time the scale factor a{t) will 



r 2c 2 + a 1 515+^ 



1+3A 
6(1+A) 



(49) 



1+3A 



The acceleration takes place when g < i.e. a(t) > [-^{A + 1)^] '^'■^+^^ 



In all the cases discussed above (i.e. for different expressions of a ), we find out 
the conditions such that g < 0. The equations (46), (48) and (49) are consistent in 
the sense that when A tends to zero both the expressions for a{tfiip) become iden- 
tical. As discussed in the end of the last section the observational constraints point 



to a tiny value of the constant A. At this small value of A the expression aitjup) 
in equation (46) is greater than that in equation (48). Since a{t) is a monotonically 
increasing function of time we get similar results from the Ray Chaudhury equation 
also in respect of the flip time which is discussed in the previous section for small 
values of A. If we consider the equation (49) for a = 1.724 as Lu's [11] choice, in 

this case A = —0.0886 (which is close to Lu's data) such that afup = ^^L^^'^^^^ ^_ 



5 Concluding Remarks 

Here we have considered the homogeneous FRW model with Modified Chaplygin 
type gas. Our analysis is based on the results of different sets of observational data. 
There is a continued debate on the exact range of the values of the exponent, a 
which generalizes the original chaplygin gas. While most observations point to the 
value of a as nearly equal to unity but existing literature abounds with examples 
of, a > 1, which incidentally may give > v^. This results in a perturbation of 
the spacetime and a perturbative analysis of the whole system shows that it favours 
structure formation. While no basic agreement is reached most workers narrow 
down the range as a is 0.9 < a < 1. Lu et al [11] gives for the MCG best fit data 
A = 0.085 and a = 1.724. In these context we have considered a = 1, a = 
and a = j^f^ which are in basic agreement with the observational analysis. Our 
findings are summarised as follows: 

1. As is well known it is very difficult to get exact form of solution of the 
field equations so we have studied graphically the variation of scale factor a{t) with 
t for different values of a. The figure shows that as a increases the rate of change 
of scale factor decreases. 

2. We have studied the key equation (10) with the help of deceleration parameter. 
From the definition of the deceleration parameter q we have calculated the fiip 
density pfup. At the larger values of a the pfup decreases, i.e., flip occurs at lower 
density or at a later time. Since the acceleration is a recent phenomena, this result 
is in agreement with the observational analysis that the value of a is nearly equal 
to the unity (a is 0.9 < a < 1). From the figure-2 it is seen that pfUp is lower for 
the higher values of a. 

3. Since our universe is accelerating our discussion emphasizes only the late stage 
of evolution. In case B we get a time dependent effective equation of state W(t). It 
gives at the late stage of evolution as a{t) — oo, W(t) — t- —1. So we asymptotically 
get p = —p from this Chaplygin type of gas, which corresponds to an empty universe 
with cosmological constant such that the equation (11) implies that the deceleration 
parameter, q tries to attain to —1. Interestingly VV(t) always remains greater than 
— 1, thus avoiding the undesirable feature of big rip. Z. K. Guo and Y. Z. Zhang [45] 
considered the new variant of CG as B{a) = Boa^"' where Bq and n are constants 
and n < 4 and Bq > 0. They finally end up with a constant equation of state 
parameter. In this case they got the EoS parameter w = —1 + n/Q, which is time 
independent. However, in our case we can avoid big rip without introducing any 
extra parameter. 

4. We have studied the velocity of sound in the Modified Chaplygin Gas model. 



Here we have discussed the possibihty of the speed of Vg is greater than the speed 
of hght. For a = 1 and A = 0, Vg is always less than Vc, but for A ^ 0, Vg exceeds Vc 
at the late stage of evolution. For a = and a = j^f^, we get Vg is always less 
than Vc- 

5. Taking first approximation of the r.h.s. of equation (10) we get the equation 
(28). For a = 1, a = and a = we get the solution of equation (28) in the 
exact form of a{t) = aosinh'^ut. We have seen that flip depends upon a. From the 
figure-7 it is seen that the flip time increases with lower value of n. Moreover the 
flip time characterized by equation (36) is found greater than that in equation (34) 
and similarly flip time for the equation (38) is greater than the equation (36). This 
finding may have some observational implications. So as a goes to unity, the higher 
value, A should vanish. This, however, is in agreement with the Fabris contention 
that recent observations point to a vanishing A. Another explanation is that if the 
value of a is greater than unity we get the negative value of A as suggested by 
Lu [11]. 

6. The whole exercise is discussed in the context of Raychaudhuri equation. As 
expected the results are in broad agreement with the previous findings. 

The main drawback of the present analysis is that we have not been able so far to 
constrain the model parameters with the help of observational data as is customary 
in relevant works in this field. It would also be a nice idea to use redshifts in place 
of cosmic time in most of the equations particularly in drawing the graphs. That 
would have been more consistent with the current nomenclature. Both the issues 
will be addressed in our future work. 
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